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Abstract 



We consider a two-way relay network, where two source nodes, SI and S2, exchange information through 
a cluster of relay nodes. The relay nodes receive the sum signal from SI and S2 in the first time slot. In the 



O 

^-j- . second time slot, each relay node multiplies its received signal by a complex coefficient and retransmits the signal 

to the two source nodes, which leads to a collaborative two-way beamforming system. By applying the principle 



of analog network coding, each receiver at SI and S2 cancels the "self-interference" in the received signal from 
the relay cluster and decodes the message. This paper studies the 2-dimensional achievable rate region for such a 
two-way relay network with collaborative beamforming. With different assumptions of channel reciprocity between 
the source-relay and relay-source channels, the achievable rate region is characterized under two setups. First, with 

> 

reciprocal channels, we investigate the achievable rate regions when the relay cluster is subject to a sum-power 

00 ' 

. constraint or individual-power constraints. We show that the optimal beamforming vectors obtained from solving 



the weighted sum inverse-SNR minimization (WSISMin) problems are sufficient to characterize the corresponding 
achievable rate region. Furthermore, we derive the closed form solutions for those optimal beamforming vectors and 
consequently propose the partially distributed algorithms to implement the optimal beamforming, where each relay 
node only needs the local channel information and one global parameter. Second, with the non-reciprocal channels, 
the achievable rate regions are also characterized for both the sum-power constraint case and the individual-power 
constraint case. Although no closed-form solutions are available under this setup, we present efficient numerical 
algorithms by solving a sequence of semi-definite programming (SDP) problems after semi-definite relaxation (SDR), 
where the optimal beamformer can be obtained under the sum-power constraint and approximate optimal solutions 
can be obtained under the individual power constraints. 
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I. Introduction 

Cooperative communication has been extensively studied in past years, where various cooperative relaying 
schemes have been proposed, such as amplify-and-forward (AF) (H, decode-and-forward (DF) O, compress-and- 
forward (CF) (3), and coded-cooperation H. Among these schemes, due to its simplicity, the AF-based relaying 
is of the most practical interest, where multi-antenna relay beamforming has also been explored to achieve higher 
spatial diversity (5J. In certain resource constrained networks, such as sensor networks, the node size is limited 
such that each node could only mount a single antenna |6). In order to exploit the multi-antenna gain in such size- 
limited cases, collaborative relay beamforming strategies have been developed where the relaying nodes cooperate 
to generate a beam towards the receiver under sum or individual power constraints EMU. 

As an extension to the AF-based one-way relaying scheme, the AF-based two-way relaying scheme (91] is based 
on the principle of analog network coding (ANC) ifTOTl to support communications in two directions. Traditionally, 
two-way relaying avoids the simultaneous transmissions of two source terminals, and requires four time-slots to 
finish one round of information exchange between them. On the contrary, the two-way relaying scheme proposed 
in iflOl allows the relay to mix the data and amplify-and-forward it, where the two terminals exploit the underlying 
self-interference structure. By doing so, the amount of required transmission time-slots is reduced from four to two 
and the overall network throughput is thus improved. There are several other works discussing such two-way relay 
systems. In particular, the authors in (9[ characterized the maximum achievable rate region for the two-way relay 
beamforming scheme by assuming a single relay node equipped with multiple antennas and two source nodes each 
equipped with a single antenna. As a counterpart of the work in (9|, the decode-and-forward two-way relaying 
has been studied in iTTTTl and the performance bounds are given in lfl2ll . The authors in (T3l studied the AF-based 
two-way relay with collaborative beamforming, where the focus is to minimize the total transmit power of the 
source nodes and the relay cluster under a given pair of signal-to-noise ratio (SNR) constraints. 

The works on characterizing the rate region of two-way relaying has also been done in (T4ll . where authors 
considered the collaborative beamforming case. However, all of existing works only obtained numerical solutions. 
The missing of closed-form solutions leads to difficulties in designing efficient algorithms due to the lack of insight 
into the structure of the optimal beamforming vectors. Thereby, in this paper, we try to seek the closed-form solutions 
for the optimal beamforming vectors to characterize the maximum achievable rate region and correspondingly 
propose efficient distributed algorithms. Our work differs from the work in O from two main aspects. First, we 
assume a cluster of single-antenna relay nodes and consider collaborative two-way relay beamforming rather than 
the multiple-antenna single-relay beamforming. Due to the distributed feature, we will study the case where each 
relay node has an individual power constraint in addition to the case where all relay nodes are subject to a sum- 
power constraint. Second, we present closed-form solutions for the optimal beamforming vectors when we have 
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reciprocal channels. For the non-reciprocal channels, the achievable rate regions are also characterized. Although 
no closed-form solutions are available under this setup, we present efficient numerical algorithms by solving a 
sequence of semi-definite programming (SDP) problems after semi-definite relaxation (SDR), where the optimal 
beamformer can be obtained under the sum-power constraint and approximate optimal solutions can be obtained 
under the individual power constraints. 

The rest of the paper is organized as follows. In Section [TTJ we introduce the system model and define the 
achievable rate region. In Sections HIT] and |W1 we characterize the achievable rate regions with two different 
assumptions on the channel reciprocity. Sub-optimal schemes with lower complexity are discussed in Section IVl 
Numerical results are presented in Section |VTJ with conclusions in Section IVIII 

Notations: We use uppercase bold letters to denote matrices and lowercase bold letters to denote vectors. The 
conjugate, transpose, and Hermitian transpose are denoted by (•)*, (-) T , and (•) , respectively. The phase of a 
complex variable a is denoted as Za. We use tr(-) and rank(-) to represent the trace and the rank of a matrix, 
respectively. A diagonal matrix with the elements of vector a as diagonal entries is denoted as diag(a). A y means 
A is positive semi-definite, a y b means ai > hi component-wise, and stands for the Hadamard (elementwise) 
multiplication. 

II. System Model 

As shown in Fig. [T] we consider a collaborative two-way relay system consisting of two source nodes SI and 
S2, each with a single antenna, and a relay cluster with K single-antenna relay nodes itVs, i = 1, ••• ,K. No 
direct links between SI and S2 exist. The forward channels from SI and S2 to relay node i axe denoted as h± t i and 
h^u respectively, while h r Xi and h r 2i denote the backward channels from relay node i to SI and S2, respectively. 
All the involved channels are assumed to take complex values and remain constant during one operation period. In 
addition, all channel state information is revealed to SI, S2, and the design/control center where the beamforming 
solution is solved. 
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The two-way relaying takes two consecutive equal-length time-slots to finish one round of communication between 
SI and S2 via the relay cluster with perfect synchronization assumed among SI, S2, and Ri, i = 1, • • • , K. In the 
first time-slot, SI and S2 transmit their signals simultaneously to the relay cluster; the i-th relay node receives the 
mixed signal ij(n), which is expressed as 

ti(n) = h^iSi(n) + h 2 ^s 2 (n) + Vi{n), (1) 

where s\(n) and S2{n) are the transmitted symbols at time index n; and Vi(n) is the receiver noise at relay node 
i, which is assumed to be circularly symmetric complex Gaussian (CSCG) with zero mean and variance af. In the 
second time-slot, upon receiving the mixed signal, relay node i multiplies a complex coefficient Wi and forwards 
the signal, which is given as Ui(n) = Witi(n). At the source node terminals, SI and S2 receive the sum signals 
from all the relay nodes, which are respectively given as 

K 

yi( n ) = y^h r i,i u i{ n ) + ^i( n )> ( 2 ) 



i=l 
K 

i=l 



2/2 (") = 2]h r 2i Ui(n) + z 2 (n), (3) 



where zi{n) and z 2 (n) are the noises at SI and S2, respectively, which are assumed to be CSCG with zero mean 
and variances a 2 si and <r| 2 , respectively. Since SI and S2 know their own transmitted signals, and s 2 (n), 

respectively, they could subtract the resulting self-interference terms Yli=i ^i i w ihi,i s l ( n ) an d Yld=i ^ 2 iWih 2t iS 2 (n) 
from the received signals, respectively. Accordingly, the remaining signals for SI and S2 are 

A 



yi( n ) = ^2 [K,i w i h 2,iS 2 (n) + hl ti WiVi(n)] + z 1 (n), (4) 

i=l 
K 

V2(n) = y~] [h^Wjh^jSiin) + h r 2i WjVi{n)\ + z 2 (n). (5) 

i=l 

Therefore, for a given w = [w±, ■ ■ ■ ,wk] T the maximum achievable rates for the end-to-end link from S2 to SI 
and from SI to S2 are respectively given as 

i?l = - log 2 1 + ; JT- , (6) 



2 S2 \ <4i + w H Aiw 



R 2 = -log a 1+ 2 2 1 ' , (7) 

where fi = hi0h^, f 2 = h 2 0h^ with h; = [h^i, ■■■ , h^x] 7 and h[ = [h\^ • • • , h r iK ] T , i = 1, 2. In addition, Ai = 
diag[|/i'J i\ 2 cr\ , • • • , \h[ k\ 2(7 k]> ^2 = diag[|/i 2 i| of, • • • , {h^ k\ 2(j2 k\> Psi an d Ps2 are the maximum transmit 
powers at SI and S2, respectively, and the factor 1/2 is due to the use of two orthogonal time-slots for relaying. 
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Accordingly, we can define the set of rate pairs achievable by all feasible beamforming vector w's as 

K= |J {(n,r 2 ) :n <Ri,r 2 <R 2 }, (8) 

where the feasible set Q w can be defined by either a sum-power constraint or individual-power constraints. 
Specifically, when the sum-power constraint is considered, we have Q w = {w : pr(w) < Pr}, where Pr is 
a scalar power limit and pr(~w) is the sum-power of the relay cluster given the beamforming vector w. When 
individual-power constraints are considered, we have Q w = {w : p_r(w) ^ Pr}, where p_r(w) is a vector of 
individual transmit powers, Pr is a vector with its elements denoting the power constraints for individual relay 
nodes, and < is element-wise. 

When time-sharing between different achievable rate pairs is considered, the achievable rate region is then defined 
as the convex hull over the set of 1Z. 

Definition 1: The achievable rate region O is the convex hull over the set of achievable rate pairs 1Z, i.e., 

= H CVK (K), (9) 

where H cvx (-) is the convex hull operation. 

The goal of this paper is to efficiently characterize the achievable rate region O. According to different assump- 
tions on the channel reciprocity between the forward and backward channels, we first study the reciprocal case, 
and then study the non-reciprocal case. 

III. Reciprocal Channel Case 

In this section, we assume that the forward channels from each source node to the relay nodes are reciprocal to 
the backward channels from the relay nodes to each corresponding source node, i.e., h± t i = h r X i and hi,% = h r 2i , 
for i = 1, • • • , K, which usually holds for a time-division-duplex (TDD) relaying system. In this case, it is obvious 
that when Zwi = —(Zhi^ + Z/i2,i), for i = 1, • • • , K, both rates given by (J6J) and (O are maximized for a given 
set of |wi|'s. Thus, we only need to further find the optimal amplitudes for the elements in w. Let Xi = \wi\ and 
fi = |^i,i|l^2,i|; we rewrite © and ©, respectively, as 

where x = [|t£>i|, • • • , |wa"|] T and f = [|/ii,i||/i2,i|> ■ ■ ■ ; |/ii,ic||^2,ftr|] • in order to obtain O, we need to characterize 
the Pareto boundary of 1Z. A common method is via solving a sequence of weighted sum-rate maximization 
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(WSRMax) problems |fT51 , each for a different non-negative weight vector (A, 1 — A), < A < 1, as follows 



max — log 2 1 + 



X 



2 &2 \ (7m + x T Aix 



51 

T 



1-A,„„ , , i%i(f x 



i2 



+— - log 2 1 + T . (12) 

2 \ a z S2 + A 2 x J 

s.t. x G (13) 

Unfortunately, we cannot derive the closed-form solution for the WSRMax problem. However, from ([Tol l and (fTTT) . 
we see that the received SNRs at SI and S2 are 

SNRl = f S f . (14) 
a z si + x 7 Aix 

SN R 2 = f^f , (15) 
cr| 2 + x-' A 2 X 

respectively, where their numerators differ by only a scalar constant. As shown later, for each given weight vector 
(//, 1 — /i), < fj, < 1, we could easily find a closed-form solution for the following weighted sum inverse-SNRs 
minimization(WSISMin) problem 

mm ^!i+x r A lX 4 2+ x^A 2 x 

^ P 52 (f T x)2 1 W P S l(f T x)2 

s.t. x G O ra . (17) 

Hence, we could quantify the Pareto boundary for the inverse-SNR region. Based on this observation, together with 
the fact that there exists a bijective mapping between an inverse-SNR pair and a rate pair, we are inspired to probe 
the question on whether we could construct the achievable rate region O from the easily obtainable inverse-SNR 
region Q. In the following, we first introduce some definitions related to the inverse-SNR region and then show that 
we indeed can construct the achievable rate region O from the inverse-SNR region, based on the set of closed-form 
solutions for a sequence of the WSISMin problems. 



A. Characterizing the Achievable Rate Region 
At first, we introduce some definitions. 

Definition 2: Consider a bijective mapping U : (x, y) \-t (| log 2 (l + 1/x), \ log 2 (l + 1/y)) with (x, y) G R++', 

1 In an independent work of [16], the authors also applied the approach of inverse-SNR. However, their motivation and objective are 
completely different from ours, where they use the sum-MSE as the objective and use the approximation MSE « 1/SNR to transform 
the sum-MSE into the sum of inverse-SNR. Hence, their results only hold for high SNR cases. In addition, they apply the beamformer that 
corresponds to the minimum sum-MSE directly to maximize the sum-rate. The optimality is suspicious since the beamformer that minimizes 
the sum-MSE is not necessary the beamformer that maximizes the sum-rate. 
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then the set of achievable inverse-SNR pairs X is defined as 

X = {(t 1 ,t 2 ) :W(ti,t 2 ) e ll}. 

For regions 1Z and X, we are particularly interested in their Pareto boundaries, which are defined as follows. 

Definition 3: The Pareto boundary of 1Z is defined as V = {{r\,r2) : (ri,r2) £ 7£, [(n>?"2) + K] = 
(^l,^)}, and the Pareto boundary of X is defined as 6 = {(^1,^2) : ^2) G X, [(ii,^) — K] X = (ii,t2)}> 
where K = R/j_ is a non-negative cone |[T5l . 

Definition 4: Define the points obtained by solving the WSISMin problem with a given weight vector (jj,, p) as 

S(n,l) = {(ti,t 2 ) : min ^1 + /2£ 2 }, (18) 

(ti,*2)ex 

where p, = I — p,. 

Definition 5: The set of points that can be obtained from a sequence of WSISMin problems is given as 

S(X)= |J SO,J). (19) 

0</x<l 

In order to show that we can construct O from 5(X), where O can be obtained by convex hulling over V, we 
need to prove two things: (1) A point in B could be mapped to a point in V and vice versa, which means U(B) = V; 
(2) The points in B that cannot be obtained by WSISMin are mapped to the points in V that are unnecessary for 
constructing O by convex hulling over V, i.e., U(B \ S(X)) C V non , where V non denotes the set of points in V 
that are unnecessary for constructing O by convex hulling over V. With the above two statements hold, it is easy 
to see that V \ V ncm C U{S(I)), i.e., the points obtained by WSISMin suffice to construct O. 

Proposition 1: The Pareto boundary of the inverse-SNR region can be mapped to the Pareto boundary of the 
rate region 1Z by mapping U as given in Definition |2j and vice versa, i.e., 

V=U(B). (20) 

Proof: See Appendix in I VIII- A I ■ 

Proposition 2: The image of B \ S(Z) is not necessary for constructing the achievable rate region O. 

Proof: See Appendix in IVIII-BI ■ 
Theorem 1: The points in S(T) are sufficient to construct the achievable rate region O. 

Proof: Since V = U(B) and U(B \ S(X)) in V is not necessary for constructing the achievable rate region O, 
it is easy to see that U{S{X)) suffices to construct the achievable rate region O given that O is obtained by convex 
hulling over V. ■ 



Since we have shown that SiT) suffices to construct the achievable rate region O, instead of studying the problem 
in ([121 and ( fT3T ), we now study the solutions of the WSISMin problems in the following. 

B. Collaborative Beamforming under Sum-power Constraint 

In this subsection, we consider the case where the relay cluster has a sum-power constraint. The total transmit 
power of the relay cluster is 

K 

PR = (\ x i h hi\ 2p Sl + \xih 2> i\ 2 P S 2 + \xi\ 2 o?) (21) 

i=l 

= x ff Dx, (22) 

where we have D = diag[|/ii 5 i| 2 P 5 i + \h 2) \\ 2 Ps2 +af,--- ,\hi, K \ 2 P,si + \h 2 ,K\ 2 Ps2 + cr 2 K }. According to the 
discussion in the last subsection, to quantify the rate region is equivalent to seeking the optimal solutions for the 
WSISMin problems given the sum-power constraint as follows: 

min )u/SNRi + /2/SNR 2 (23) 

X 

s.t. x T Dx < P R , (24) 

where < [i < 1 is the weight. First of all, the optimal x* must satisfy x* T Dx* = Pr; otherwise, we can always 
scale up x such that the objective function is decreased. When fi = or 1, this problem degrades to find the optimal 
beamforming vector for collaborative one-way relay, which has been extensively studied in Q, Q, and ifToll . 

Given the SNRs from (0 and (Q3]) and the fact of x* T Dx* = Pr, the optimal x for (|23]) should be the solution 
of the following problem: 

. x T [z.D/P fl + ^/P S2 A 1 + /i/P5iA 2 ]x 

mm — , (25) 

x T ff r x 

where v = ^o- 2 sl /Ps2 + P-c^/Psi- The above problem is equivalent to 

x r ff T x 

m f C x^[^D/P ii + M /P S2 A 1 + /i/P5iA 2 ]x' ( } 

which can be written in the form of Rayleigh-Ritz ratio ifTTl and the optimal solution is thus given as 

x* = er^f/iir^fn, (27) 
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^ 2 +Psi|/ll,i| 2 +Ps2|/l2,. 
2 +P S l|/il.,| 2 +P S2 |/l2, l | 2 



2 



SNRi — — — , (33) 

°S1 Z^i=l CT 2 



p 5i fe£i I^Pi^A/^+p^i^ ^+p^i^ai 
SNR 2 = ^ V ,,. ' , — '—. (34) 

°S2 "T Z^i=l <r 2 +P sl |/i 1 , 1 | 2 +P S2 |h 2 , I | 2 



where 



T = diag 



(28) 



A = ^ 2 + P5iM 2 + ^2|/i2,i| 2 , (29) 
?7i = of ( |/il,i| 2 /U/P S l + |/i2,j| 2 £/-Ps2 ) , (30) 



and £ is a scalar such that x* T Dx* = Pr. By searching over all /i's, we derived a set of x*'s and hence we could 
compute a set of rate pairs by injecting (1271) into (fTOl) and (fTTb . The achievable rate region is then obtained by 
convex-hulling over such a set of rate pairs. 

Partially distributed implementation: After exhausting all /z's, we set up a lookup table for rate -pairs indexed 
by fj,. During normal operations, the control center first looks up the table to decides the appropriate (jl such that 
SI and S2 achieve a desirable rate pair; and it broadcasts \i and the global constant £/||r -1 f||, while Pgi, Ps2, 
o~si, &S2 are constant and assumed to be known at all the relays. Upon receiving the broadcast message from the 
control center, each relay node determines the optimal Wi from its local information hij and h,2,i, which is given 
as 

e IMIM xth^+zh^) (31) 

||r- 1 f||^/P R + 7? l 

C. Collaborative Beamforming under Individual-Power Constraints 

In the previous subsection, we assume that the relay cluster has a sum-power constraint. In practice, each relay 
may have its own power constraint due to the individual power supplies. The transmit power at relay i is given as 

PR,i = \vi\ 2 (\hi,i\ 2 P S i + \h 2 ,i\ 2 Ps2 + erf), (32) 

where pn^ < pi, with pi is the maximum allowable power for relay node i. Equivalently, we could set pn^ = afpi 
with < ai < 1 as a new design variable. Correspondingly, the received SNRs can be rewritten as (l33l) and (l34l) 
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Let 



Hi 
H 2 

9i 



diagfo-ipil/i^il 2 , • • • ,cr 2 K p K \hi :K \ 2 }, 
diag[<TiPi|/i 2i i| 2 , • • • ,cr 2 K p K \h2,K\ 2 ]i 



/Pi\hi,i\\h2,i\/ y/Qi 



(35) 
(36) 
(37) 



We can recast (l33l and (l34l as 



SNRi 
SNR 2 



-Pg2Q gg Q 
a 2 +a r HiD" 1 a : 



PgiQ i gg^ a 
<r| 2 + a T H 2 D- 1 a i 



(38) 
(39) 



respectively, where ■< a ^ 1. The WSISMin problem for the individual-power constraint case is now given as: 

v + a T (HiD~VAPs2 + HaD-^/PsOa 



mm 



which is equivalent to solve 



max 



a 1 gg 1 a. 



T T 

a gg a 



(40) 



o^i v + a T (H.iY)- 1 fi/Ps2 + HaD-i/i/PsOa ' 



(41) 



For notation simplicity, let 



* = [(Hi/VP S 2 + H 2 /l/Psi)D-7 
g = g/v^j 



- 1 VI 1/2 



(42) 
(43) 



where * is diagonal with its diagonal elements denoted as ipi, i = 1, • • • ,-ftT. Then the above problem becomes 

<g,<*> 2 



max 



o^a^i 1 + I *a| 



(44) 



For each given /i, (1441 can be solved analytically by following the results in Qjjj. Before we present the solution, 
we first define fa = cn/ipf for i = 1, ■ ■ ■ , K and <^k-+i = 0. Then we sort fa as <f> Tl > fa 2 > • • • > TK > 
Moreover, let = 1 tS"' =1 -^ Tm and define the j-th element of the vector a.^ as 

1, j = Ti,---r k 

^k4>j, 3 = Tk+l-, ■ ■ ■ T K 

Then the solution for (1441) is given by following theorem. 



(fc) 



(45) 



Theorem 2: The solution of (l44l) is given by (|45T ). where fc* is the smallest fc such that Afc < <j> T ^ +1 - 

Proof: This result directly follows the results in (H. 
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Partially distributed implementation: Besides the value of fi, the control center only needs to broadcast A&* at 
each operation period. Each relay node then determines <f>i with its local information. If (f)^ 1 < A&*, the relay 
node transmits at its maximum power. Otherwise, it transmits with power {Xk-^ifpi, i- e -> the optimal Wi = 
a [ k \/pie~^ zhlA+zh2 ''\ where a[ k is given in (|45T ). From the solutions, we see that in general some relay nodes 
may not transmit with maximum transmit power. 

IV. Non-reciprocal Channel Case 

In the last section, we have discussed the case where the uplink and downlink channels are reciprocal. In this 
section, we discuss the case where the uplink and downlink channels are non-reciprocal, which may be the result 
of deploying frequency-division-duplex (FDD) system. 

Due to the lack of channel reciprocity, the approach taken in the last section does not apply here. In order 
to characterize the boundary of the region 1Z, as we discussed before a commonly used method is to solve the 
following 

A / P 52 |f 2 T w| 2 
max — log 2 I 1 + 



w 



2 V <?% +w ff Aiw 



si 



2 V ffg 2 + w JJ A2W/ 

s.t. w G Q w , (47) 

for each given weight vector (A, 1 — A). However, the above problem is non-convex since the objective function is 
not a concave function. To efficiently quantify the rate region, here we resort to an alternative method called the 
rate-profile method |9), formulated as 

max R sum (48) 

s.t. - log 2 1 + 2 7 1 2 H ' > KR sum , (49) 

V °si + w Aiw/ 

\ log 2 (l + ^^h! ) > *Rsu m , (50) 

w£n w , (51) 
where R sum is the sum rate given a rate profile vector [«, «] with < k < 1 and R = 1 — k. Let Fi = f*f^, 
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F2 = f|fj, and X = ww^. The above problem is equivalent to 



max R sum (52) 

1 / P52tr(F 2 X) \ 

2 v 4i + tr ( A i x )/ 

1 i / 1 j Psitr(FiX) 
2° g2 V 1+ a| 2 + tr(A 2 X) 

X G Ox, (55) 



- log 2 ( 1 + — — ; ^ K ^ ) > KR sum , (54) 



X ^ 0, (56) 
rank(X) = 1, (57) 

where the last constraint rank(X) = 1 comes from the fact X = ww^, and Ox = {X : X = ww^,w £ O w } 
and Q w is defined after ([8]). According to different assumptions on the power constraint, the above problem can be 
further converted into different semi-definite programming (SDP) problems after semi-definite relaxation (SDR). 

A. Sum-power Constrained Case 

In this subsection, we assume that the relay cluster operates under a sum-power constraint Pr. Given the sum- 
power constraint, the power constraint in (1551 ) can be replaced by tr(DX) < Pr, where D = diag[|/ii i| 2 Pgi + 
|^2,i| 2 Ps2 + • • • , \h\,K\ 2 Psi + \h2,x\ 2 Ps2 + c 2 ^]- Since the rank-one constraint is not convex, the problem is 
still not a convex problem and hence may not be efficiently solvable. To address this issue, let us first remove the 
rank-one constraint and consider the following relay power minimization problem for given set of k and R sum = r: 



min 


tr(DX) 


X 




s.t. 












x y o, 



(58) 



P S2 tr(F 2 X) 



1 + V h > itr, (59) 
4 1 +tr(A 1 X)' 

' Pgitr(FiX) 
a 2 S2 + tr(A 2 X) 

(61) 



which is equivalent to 



min tr(DX) (62) 



P S2 tr(F 2 X) 

s-t 2 | + { \ "v" \ ^ 7i» (63) 
a z si +tr(AiX) 

Psitr(FiX) 

> 72, (64) 



<>S2 + tr(A 2 X) 

X h 0, (65) 
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where 71 = 2 2ra " — 1, 72 = 2 2Kr — 1, and they can be considered as the SNR constraints for SI and S2, respectively. 
Since a 2 si + tr(AiX) > and cr| 2 + tr(A 2 X) > 0, we could rewrite the above problem as following SDP problem: 

min tr(DX) (66) 

s.t. tr[(P 52 F 2 - 7iAi)X] > 71^!, (67) 

tr[(P 5 iFi - 72 A 2 )X] > 72( t| 2 , (68) 

X >: 0. (69) 

Denote the optimal value of the above problem as p* R , which is the minimum sum-power required by the relay 
cluster to support the target SNRs 71 and 7 2 for SI and S2, respectively. If p* R < Pr, then (71,72) must be an 
achievable SNR pair. Otherwise, 71 and 7 2 are not achievable. Based on this observation, we propose the following 
bi-section algorithm such that the problem (T52l) without rank-one constraint can be solved by solving a sequence of 
convex power feasibility problems, with the assumption that we know an upper bound for R sum , denoted as r max . 
Algorithm 1: 

. Initialize r low = 0, r up = r max . 

• Repeat 

1) Set r <- ±(r Zoi0 + 

' up)- 

2) Solve problem (I66l)-(l69l with the given r. 

3) Update r with the bi-section method fi31l : If p* R < Pr, set ri ow = r; otherwise, r up = r. 

• Until r up — ri ow < e, where e is a small positive accuracy parameter. 

The rate upper bound r max can be derived as follows. We first decouple the two-way relay channel into two 
one-way relay channels and obtain a rate for each one-way relay channel. Denote the larger rate as f. Then r max 
can be set as 2f. The one-way collaborative relay beamforming with sum-power constraint is well-studied, and the 
rate can be derived from the results in 0. 

1 ) Rank-one solution: The resulting optimal solution X opt obtained from Algorithm 1 may not be of rank-one 
due to the SDP relaxation, which means that X opt may not lead to an optimal beamforming vector w. However, 
since there are only two linear constraints (f6Tb and d68l) , it has been shown in and |[T8l that an exact rank-one 
optimal solution can always be constructed from a non-rank-one optimal solution. The transformation techniques 
developed in [9] and |fT8l can be used to obtain the rank-one solution. Note that the beamforming solution for the 
non-reciprocal channel case is fully centralized, which cannot be implemented in a partially distributed fashion. 
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B. Individual-Power Constrained Case 

In the previous subsection, we have discussed the sum-power constrained case where the non-convex rate 
maximization problem is converted into a sequence of convex sum-power minimization problems. In this subsection, 
we put a stricter limitation on the relay power by assuming that each node has its individual power constraint. 
In this case, following a similar SDR technique to that in the previous subsection, the optimization problem with 
individual power constraints can be cast as 

max R sum (70) 
P S2 tr(F 2 X) 

S-t. 2 , w A y n > 71, ( 71 ) 

a z si + tr(AiX) 
°S2 + tr(A 2 X) 

D M X M <P^ i = 1, ■ ■ (73) 
X >: 0, (74) 

where Dj j and Xj j are the i-th diagonal elements of D and X, respectively. The transmit power at node i amounts 
to Dj /Xj j and the individual power limit at node i is P % R . However, we cannot translate the above problem into 
a sequence of power feasibility problems as given in the last subsection, since we now have K individual power 
constraints rather than a single sum-power constraint for the whole relay cluster. Alternatively, we aim at solving 
a sequence of the following problem via bi-section search over r. 

max r (75) 

X,r 

s.t. tr[(P S2 F 2 - 7 iAi)X] > 7l< 7| l5 (76) 

trKPsiFi - 72 A 2 )X] > l2 ci 2 S2 , (77) 

X(m) < P l R /T>(i,i), i = l,-- - ,K, (78) 

X h 0. (79) 
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The above problem is convex over X at each given value of r. Let r* be the maximum value obtained by solving 
(|75T ). For a given value of r, we solve the following feasibility problem 

Find X (80) 

s.t. tr[(P 52 F 2 - 7iAi)X] > tktI!, (81) 

tr[(P S iFi - 72 A 2 )X] > 7 2 4 2 , (82) 

X(i,i) <Py-D(i,i), i = l,---,K, (83) 

x y 0. (84) 

If it is feasible, we have r < r* and the corresponding rate is achievable. Otherwise, we have r > r* and the 
corresponding rate is not achievable. Based on this observation, we apply bi-section search over r to solve the 
problem in (|75T ). where we solve a convex feasibility problem of (l80l at each step. We start with an interval 
[0, r max ] that contains the optimal value r* where r max can be obtained in a similar way as that for the sum-power 
constrained case, and run the following algorithm. 
Algorithm 2 

. Initialize r iow -0, r up = r max . 

• Repeat 

1) Set r <- \{r low + 

1 up)- 

2) Solve the feasibility problem given by (f80l-(l83l with given r. 

3) Update r: If the problem is feasible, set ri ow = r; otherwise, r up = r. 
. Until r up - ri ow < e. Then r* = ri ow . 

1 ) Rank-one solution based on randomization: Similar to the sum-power constrained case, the solution of X at 
the end of Algorithm 2, denoted as ~K opt , may not be rank-one. However, since there are K + 2 linear constraints 
here, we cannot apply the rank-one decomposition technique in lfT8l . which require the number of linear constraints 
to be less than or equal to 3. Fortunately, various techniques have been developed ||T9l to generate good rank- 
one approximate solutions to the original problemy. One such efficient approach is based on randomization |[T9l : 
using X opt to randomly generate a set of candidate weight vectors, {w^}, from which the "best" solution for the 
beamforming vector w is selected. There are three ways of generating {w^} as presented in |[T9l . In order to satisfy 
the individual power constraint, we adopt the routine named randB in |[T9l . Specially, let be the vector whose 
elements are independent random variables uniformly distributed on the unit circle in the complex plane, i.e., its 
z-th element [ej]j = e 7 where 6i/s are independent and uniformly distributed over [0, 2tt). We choose w; such 

2 The randomization technique only provides approximate solutions. Hence, the corresponding rate region is not exact. 
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that its z-th element [wj]; = y/[X.opt]u [ e l]i- As we see, |[w/]j| 2 = [X opt ]jj; hence the individual power constraint 
can be satisfied. 

For each x'^wiwf, we associate each with a value v(wi), 

v(wi) = max ( 1 - tr[( jPsi ^ 1 - ^-)w ; wf ], 

l-tr[(^-4l)w,wf]Y (85) 
71^51 a Si / 

which reflects how much the constraints are violated. The "best" weight vector among the candidate vectors is the 
one that has the minimum v(wi), i.e., 

I* = arg min v(wi), (86) 
w* = w>. (87) 

V. Sub-optimal Schemes 

In this section, we propose some suboptimal schemes with lower complexity for implementation than the optimal 
ones established in the previous sections. 

A. Reciprocal Channel Case 

In the reciprocal channel case, at first the transmit phases #j's at the relays are matched to the channels as 
Oi = — (Z/ix j + Z/i 2> j). Then with the sum-power constraint, we propose the sub-optimal equal power beamforming 
scheme where each relay transmits with equal power. With the individual-power constraints, we propose the max- 
power beamforming scheme where each relay transmits with its maximum power. 

1) Equal-power beamforming: All the K relay nodes transmit with the same power Pr/K; 9i's and Xj's for 
i = 1, • • • , K, are given as: 

9i = -(Zh hi + Zh 2 ,i), (88) 

Xi = ^ KiPs^tf + Pstlhztf + a?)- (89) 

2) Max-power beamforming: Each relay transmits with its maximum allowable power Pr^\ 9i's and Xj's for 
i = 1, ■ ■ ■ ,K, are given as: 

Oi = -(Zh lii + Zh 2 , i ), (90) 



Pr, 

Psi\hi,i\ 2 + Ps2\h2,i\ 2 + of 



x i ~ \l r, ir_ 19 , n it. ,|2 , _2' (91) 
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These sub-optimal schemes enjoy implementation simplicity since each relay only requires the local channel 
information hij and ^2,1 to decide the transmit phase and Xj. 

B. Non-reciprocal Channel Case 

For the non-reciprocal channel case, since the transmit phase cannot be matched to the two-directional channels 
simultaneously, we propose a sub-optimal scheme that greedily chooses the transmit phases. Specifically, each relay 
chooses the transmit phase to be either Z/ii^ + /-h r 2 i or Z/12,?. + ^-h\ p whichever maximizes its own contribution 
to the overall SNRs at SI or S2 without considering any other relays' contributions, i.e., we pick one of the above 
two phases that maximizes the following quantity: 



where x% is the transmit amplitude. To determine Xj's, we adopt equal-power beamforming for the sum-power 
constraint case and max-power beamforming for the individual-power constraint case, which are given in (l89l and 
(I9TT) . respectively. 



In the section, we present numerical results to quantify the achievable rate region for the two-way relay network 
with collaborative beamforming. We assume that the relay cluster consists of 5 nodes; the channel coefficients h\i 
and h\ ^ i = 1, ■ ■ ■ , K, are independent CSCG variables with distribution CAf(0, 1); the channel coefficients h^i 



are assumed to have unit variance in the simulations. We change fj, from to 1 with step 0.1 and obtain 11 Pareto 
boundary points. For each point, we run 100 channel realizations to measure the expected performance. We then 
do convex hulling over these points. 

First, we investigate the achievable rate region. Both reciprocal channel and non-reciprocal channel cases are 
discussed under a sum-power constraint and individual-power constraints, respectively. When the channel is recip- 
rocal, we set hi t i = h r Xi and = hT 2 4 for i = 1, ■ ■ ■ , K. As shown in Fig. 12 the solid curves represent the 
reciprocal channel cases, with the outer one denoting the sum-power constraint case and the inner one denoting 
the individual-power constraints case; the dashed curves represent the non-reciprocal channel cases, with the outer 
one denoting the sum-power constraint case and the inner one denoting the individual-power constraints case. For 
the sum-power constraint case, the relay power Pr = 10 W while the transmit powers Pgi = Ps2 = 1 W. For the 
individual-power constraints case, the relay power constraints are given as 2.5, 3, 0.5, 1, 3 W (noises are assumed 
to have unit power in Watt), which is summed up to 10 W. We use CVX, a Matlab-based optimization software 




(92) 



VI. Numerical Results 



and h T 2i , i = 1 



, • • • , K, are also independent and distributed as CAf(0, 1). The noises at the relays and source nodes 
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Fig. 2: Achievable rate regions for reciprocal and non-reciprocal channels with sum-power constraint and individual-power 
constraint, respectively. Transmitter powers: Psi=Ps2=i W, relay network power Pr=10 W (sum-power constraint), Pr= 

[2.5, 3, 0.5, 1, 3] W (individual power constraint). 




0.2 0.4 0.6 0.8 1 1.2 1.4 

t i b/s/Hz 

Fig. 3: Achievable rate regions for reciprocal channel case under a sum-power constraint, network total power are 0, 10, 100 

W, using equal-power beamforming as sub-optimal scheme. 



|[20l . to solve the SDP problems. As we see in Fig. [2j due to the symmetry of the transmit powers and channel 
statistics, the achievable rate region O is symmetric. When P$2 = 0, the rate pairs collapse to the segment on the 
horizontal axis, which corresponds to the achievable rate for a one-way relay network where only SI transmits. 
Moreover, the rate region for the individual-power constraint case is smaller than that for the sum-power constraint 
case. This is quite intuitive since the individual-power constraint is stricter than the sum-power constraint. 

In Fig. |2j we also compare the rate regions for the reciprocal and non-reciprocal channel cases under the same 
power constraint assumption. As we can see, the maximum rate for SI in the reciprocal channel case is the same 
as the one in the non-reciprocal channel case. This is because such a maximum rate is obtained by optimizing the 
one-way link from SI to S2 without considering the link from S2 to SI. Since the one-way link from SI to S2 
consists of hi and h^, whether = hi or not does not affect the statistics of the one-way link from SI to S2. The 
same argument holds for maximum rate at S2. We also observe that the rate region for the reciprocal channel case 
is larger than that in the non-reciprocal channel case given the same settings of powers and noises. The reason is 
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Fig. 4: Achievable rate regions for reciprocal channel case under individual-power constraints, network total power are 
0, 10, 100 W, using max-power beamforming as sub-optimal scheme. 



100 w 




Fig. 5: Achievable rate regions for non-reciprocal channel case under a sum-power constraint, network total power are 
0, 10, 100 W, using equal-power beamforming as sub-optimal scheme. 

that we can match the beamforming phase to the overall channel phase (i.e., Zwi = Zh± t i + Zh^i, i = 1, • • • , K.) 
in the reciprocal channel case, while we are not able to do so in the non-reciprocal channel case. Therefore, TDD 
based system is more favorable in terms of the achievable rate region if the channel coherence time is larger than 
one operation period and the transmit-receive chain calibration ||2T1 can be properly done. Besides the rate region, 
the amount of information needs to be broadcast by the control center is significantly different. In the reciprocal 
channel case, the control center only needs to broadcast one scalar at each time slot. However, in the non-reciprocal 
channel case, the control center needs to broadcast the beamforming vector, which is a complex vector of dimension 
K. 

Second, we investigate the performance of the sub-optimal schemes in relative to the maximum achievable rate 
regions. As we see in Fig. [3] for reciprocal channels case under the sum-power constraint, the rate pairs achieved by 
the equal-power beamforming scheme, denoted as single points, are strictly sub-optimal. On the contrary, as shown 
in Fig. |4] for the individual-power constraints case, the rate pair achieved by max-power beamforming gets closer to 
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O.B 



Fig. 6: Achievable rate regions for non-reciprocal channel case under individual-power constraints, network total power are 
0, 10, 100 W, using maximum-power beamforming as sub-optimal scheme. 

the boundary when the power budget is reduced[] In Fig. [5] and Fig. [6l we consider the non-reciprocal channels and 
show the performance of equal-power beamforming and max-power beamforming with greedy phase selection as 
given in (l92b - The performance of both equal-power beamforming and max-power beamforming schemes degrades 
as Pr increases. Thereby, the sub-optimal schemes for the non-reciprocal channel case works well only when Pr 
is small. 

VII. Conclusion 

In this paper, we considered the two-way relay networks with collaborative beamforming and investigated the 

achievable rate region, which is defined as the convex hull of all achievable rate pairs. We studied both the 

reciprocal and non-reciprocal channel cases. In the reciprocal channel case, we characterized the rate region when 

the relay cluster is subject to either a sum-power constraint or individual-power constraints, respectively. It was 

shown that we could characterize the whole achievable rate region via the Pareto-optimal beamforming vectors 

obtained from solving a sequence of WSISMin problems. Furthermore, we derived the closed-form solutions for 

those optimal beamforming vectors and consequently proposed partially distributed algorithms to implement the 

optimal beamforming, where each relay node only needs its own local channel information and one global scalar 

sent from the control center. For the non-reciprocal channel case, we used the rate -profile approach to compute the 

Pareto-optimal beamforming vectors. When the relay cluster is subject to a sum-power constraint, we computed 

the optimal beamforming vector via solving a sequence of relaxed SDP power minimization problems followed by 

a special rank-one reconstruction. When the relay cluster is subject to individual-power constraints, we solved a 

sequence of relaxed SDP feasibility problems and the rank-one solution is obtained by randomization techniques. 

From the numerical results, we found that the achievable rate region is larger in the reciprocal channel case than 

3 We set the individual powers P_r = [2.5, 3, 0.5, 1, 3] W with total power equal to 10 W. When total power is changed to 1 W and 100 
W, we scale the vector proportionally. 
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dp 2 (y) 
dy 2 



dp'(y)/dx 
dy/dx 
x + 



In 2 



q (x)(x + x 2 
q(x) + q{x)' 



1 + 2x)(q(x) + ^(x) 2 ) — (q (x) + 2q(x)q (x))(x + x 1 
q(x) + q(x) 2 



> 0. 
(93) 



that in the non-reciprocal channel case. Hence, TDD-based relaying scheme is more favorable for the two-way 
relay network with collaborative beamforming. 

VIII. Appendices 

A. Proof of Proposition [7] 

Proof: We will show this by contradiction. Assume (a, b) G B but U{{a,b)) V. Then we can find another 
point (c,d) G 1Z such that c > l/21og 2 (l + 1/a) and d > l/21og 2 (l + 1/^)- According to the definition of Z, the 
point ( 2 2c 1 _ 1 , nzi-i ) ^ Thus, there exists a point in X such that 22c 1 _ 1 < a and 22rf 1 _ 1 < 6, which contradicts 
the assumption that (a,b) is a Pareto optimal point. Hence U{B) C "P. The converse that U{B) 5 can also be 
proven in the similar way. Therefore, V = U{B). ■ 

B. Proof of Proposition \2\ 

In order to prove Proposition [2j we first introduce the following two lemmas. 

Lemma 1: Suppose q(x) is a positive, decreasing, and linear function with x > 0. The bijective mapping U maps 
(x,q(x)) to {y,p(y)); then p(y) is a non-negative, decreasing, and convex function. 

Proof: Let y = log 2 (l + 1/x) and p(y) = log 2 (l + l/q(x)) be an implicit function of y, where x > 0. 
Since q(x) is positive, decreasing, and linear, we have q(x) > 0, q'(x) < 0, q (x) = 0, and hence p(y) > 0. The 
first-order derivative of p(y) is 

^ = dp{y)/dx 
dy/dx 

\(x) + q(x) 

The second-order derivative is given by (l93l . which is positive 
Thus, p(y) is a convex function of y. 



'(x) — ; _^ N2 < 0. 



According to the above lemma, the line segment Q1Q2 in Fig. |7(a) is mapped to a convex curve Q1Q2 in Fig- 



|7(b)| by U. In addition, it is easy to see that Q1Q2 + Q1Q2 ~ K, i.e., any point above Q1Q2 (for example, 
P in Fig. |7(a) 1 will be mapped to be a point below Q^Q^ (i.e., P in Fig. |7(b)[ ). 



22 



A 



(a) A straight line in X. (b) The image of in 

n. 

Fig. 7: Illustration of Lemma Q] Mapping a straight line in I to a convex curve in 7Z. 



Lemma 2: Let a point (91,(72) € bd(X) \ <S(//,X), where bd(X) denotes the boundary of region X. If 91 = 
Aii + Asi, where (t\, t 2 ), (si, s 2 ) € S(/i,X) and (ti, £ 2 ) # (si,s 2 ), we have g 2 > Ai 2 + As 2 , i.e., the point (gi,g 2 ) 
is above the line segment connecting (ii,£ 2 ) and (si,s 2 ). 

Proof: We show this by contradiction. Suppose <S(/i,X) has more than one elements for a given /i, such that 
(ti, t 2 ), (si, s 2 ) G S(^,l), and (ii,i 2 ) / Osi,s 2 ). According to the definition of 5(//,X) given by (fT8l ), we have 
/iti + /ii 2 = fis\ + /2s 2 = ™, where m is the minimum value of the weighted sum for a given fj, over all points in 
X. If q\ = \t\ + Asi and q 2 < A£ 2 + As 2 , we have 



M9l + < ^{\ti + \S\) + n(\t 2 + \S 2 ) 

= X(fiti + /it 2 ) + A(^si + /us 2 ) 



m. 



(94) 
(95) 
(96) 



If iiqi+Jlq 2 < m, it contradicts that m is the minimum value of the weighted sum for the given /i; If ^q\+Jiq 2 = m, 
it contradicts that (91, 9 2 ) is not in S([i,Z). Therefore, the lemma holds. ■ 

According to the above lemma, for a given fi, if <S(/z,X) has more than one elements, the set of boundary points 
{(<7i;<72) : (^1,92) £ bd(X) \5(;U,X),si < 91 < £1, (si, s 2 ), (£1, £ 2 ) £ <S(^,X)} must be above the line segment 
connecting (si, s 2 ) and (ti,t 2 ); and hence are not attainable by solving WSISMin. This is true for all /i's; hence if 
a boundary point is not attainable by solving WSISMin, it must be above a line segment connecting two particular 
points in S{p.,Z) for some p,. With the above two lemmas, we are ready to prove Proposition [2] as follows. 

Proof of Proposition |2j 



Proof: First we define 



A = {/i : S(/j,,I) has more than one elements}, 



(97) 



and let be the line segment (e.g., AB in Fig. |8(a)| ) with two end points from S(fi,I) for fi £ A (e.g., points A 
and B in 8(a) i. According to Lemma |2l the boundary points that are not attainable by solving WSISMin, denoted as 
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(a) Inverse-SNR region (b) Rate region 

Fig. 8: Inverse-SNR region and corresponding rate and region. 

bd(Z)\<S(Z) (here referring to curve AB in Fig. |8(a)| ), must be above l^'s, i.e., bd(Z)\<S(Z) C U^eA^M+K); an( * **■ 
follows that W(bd(X)\5(X)) C W(|J M6A (^+K)). According to Lemmaffl ^(U AteA (V+ K )) Q IUa(^(*m)- K )' 
where is a convex curve (e.g., here refers to the dashed convex curve A'B' in Fig. |8(b)[ ). Let 1^ be a line 
segment (i.e., the dot-dashed line segment A'B' in Fig. |8(b)| ) that connects the two end points of the convex curve 
lUlft). Due to the convexity of we have U{1^) -KC^-K and hence W(bd(J) \ S(X)) C U^eA^A* ~ K )- 

Notice U^£a(^+K) is sufficient for constructing O by convex hulling. Therefore, W(bd(X)\iS(X)) or bd(X)\«S(X) 
is not necessary for constructing 0. Since B C bd(X), the set B\S(I) is also not necessary for constructing C ■ 
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